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" Abstract 



Context. A new class of super-luminous transients has recently been identified. These objects reach absolute luminosities of 
M u < — 21, lack hydrogen in their spectra, and are exclusively discovered by non-targeted surveys because they are associated 
with very faint galaxies. 

Aims. We aim to contribute to a better understanding of these objects by studying SN 2006oz, a newly-recognized member of this 
class. 

Methods. We present multi-color light curves of SN 2006oz from the SDSS-II SN Survey that cover its rise time, as well as an optical 
spectrum that shows that the explosion occurred at z ~ 0.376. We fitted black-body functions to estimate the temperature and radius 
evolution of the photosphere and used the parametrized code S YNOW to model the spectrum. We constructed a bolometric light curve 
and compared it with explosion models. In addition, we conducted a deep search for the host galaxy with the 10 m GTC telescope. 
Results. The very early light curves show a dip in the g- and r-bands and a possible initial cooling phase in the u-band before rising 
to maximum light. The bolometric light curve shows a precursor plateau with a duration of 6-10 days in the rest-frame. A lower limit 
of M u < —21.5 can be placed on the absolute peak luminosity of the SN, while the rise time is constrained to be at least 29 days. 
During our observations, the emitting sphere doubled its radius to ~ 2 x 10 5 cm, while the temperature remained hot at ~ 15000 K. 
As for other similar SNe, the spectrum is best modeled with elements including O II and Mg II, while we tentatively suggest that Fe III 
might be present. The host galaxy is detected in gri with 25.74 ± 0.19, 24.43 ± 0.06, and 24.14 ± 0.12, respectively. It is a faint 
dwarf galaxy with M g — —16.9. 

Conclusions. We suggest that the precursor plateau might be related to a recombination wave in a circumstellar medium (CSM) and 
discuss whether this is a common property of all similar explosions. The subsequent rise can be equally well described by input from 
a magnetar or by ejecta-CSM interaction, but the models are not well constrained owing to the lack of post-maximum observations, 
and CSM interaction has difficulties accounting for the precursor plateau self-consistently. Radioactive decay is less likely to be the 
mechanism that powers the luminosity. The host is a moderately young and star-forming, but not a starburst, galaxy. 

Key words, supernovae: general, supernovae: individual: SN 2006oz 
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1. Introduction 

Historically, supernovae (SNe) have usually been discovered by 
monitoring bright, nearby galaxies. With a few exceptions, the 
supernovae discovered in this targeted way fit well withi n the tra- 
ditional SN classification scheme (e.g. lFilippenkd 1997b . During 
the last few years, however, rollin g searches, such as the Texas 
SN Search (lOuimbv et all l2005h an d the Catalina Real-Time 
Transient Survey (Drake " al.ll2009l) have changed our view of 
stellar explosions. In particular, it has been shown that interest- 
ing transients were previously missed exactly because they occur 
in environments different from those probed by the traditional 
SN searches. Remarkably, these include the brightest SNe ever 
recorded (hereafter super-luminous supernovae, SLSNe) with 
luminosities exceeding those of SNe la by 10-100 times. 

Many of these SLSNe are associated with SNe Iln and 
have high luminosities attributed (at least partly) t o the inter- 
action of the ejecta with a de nse H-rich CSM (e.g. Ofe k et all 
l2007t ISmith et al.ll2007l 12008b . SN 2005ap, in contrast, showed 
only weak evidence for hydrog en and yet reached a n absolute 
magnitude of Mr < -22.5 dOuimbv et alJ 12007b . A seem- 
ingly unrelated object of peculiar nature was disc overed by the 
Super nova Cosmology Project (SCP) with HST dBarbarv et al.l 
2009). SCP06F6 had unprecedented spectra and light curves, but 
the lack of a robust redshift estimate left this study inconclu- 
sive with respect to its nature (see also Chatzopo ulos et ail2 009; 

Idansicke et al.ll2009b ISoker et aUl2010h . " 

Significant progress was made when lOuimbv et al.1 d2011l) 
showed that four objects detect e d by the Palomar Transient 



Factory (PTF; iRau et all 120091: lLaw et al.1 [2009b could be 
grouped together with SN 2005ap and SCP06F6 to form a 
distinct class of H-poor SLSNe. The redshifts of these ob- 
jects were identified by the detection of the Mg II A2800 dou- 
blet and SCP06F6 was shown to be at a much higher redshift 
(z = 1.189) than originally estimated. All objects had simi- 
lar spectra, blue colors, and relatively symmetric light curves. 
Their typical absolute magnitude is M ~ —21.5 and their 
explosion mechanism remains a my stery. Possible sugge stions 
include pulsational pair-instability ( Wooslev et al. 2007 ]), the 
powering of th e ejecta by a magnetar (Kasen & Bildsten 2010; 
WooslevlboiOh. or interaction with an CSM (Chevalier & Irwin 



Wooslev ZU1U), or interaction with an LaM (Chevalier & Irwin 
201 U iBlinnikov & Sorokinall2010t IMoriva & To minaga 201 lh. 



A third class of SLS Ne, also H-poor, but different from those in 
lOuimbv et all d201 lh. might be represente d by SN 2007bi, which 
has been proposed dGal-Yame t al. 200% to be a pair-instability 
event, although this explanation is not unique ( Young et al. 2010; 



IMoriva et al.ll2010h . 

One of the objects studied by lOuimby et al.l d20 1 lh. 
SN 2 010gx, was also extensively followed by iPastorello et al] 
d2010l) . By obtaining spectra at later phases, they managed to 
demonstrate that this SLSN transitioned to a SN Ic, showing 
that there is a possible link between these energetic explosions. 
Recently, two more SNe belonging to thi s intriguing class wer e 
discovered by Pan-STARRS 1 at z ~ 0.90 dChomiuk et al.l201 lh . 
These objects showed no signs of deceleration in their expansion 
velocities during observations obtained over a period of about 
three weeks around maximum light. 

A common characteristic for H-poor SLSNe is that they are 
systematically found in faint galaxies. Indeed, to date, only three 
of these hosts have been detected; there are only upper limits on 
the others (typically > —18), suggesting that they are prob- 
ably metal-poor dNeill et al. 2011). It has been speculated that 
low metallicity might be an indispensable ingredient to produce 
SLSNe. 



One reason that our understanding of these objects is lim- 
ited is that the available data are sparse and incomplete: obser- 
vations are often obtained in 1 or 2 neighboring filters, and when 
multi-color light curves are available, they cover only part of the 
SN evolution. The study of SLSNe is a relatively new topic and 
any complementary dataset (especially covering the critical early 
phases) constitutes a valuable contribution to the field. 

SN 2006oz was discovered by the SDSS-II SN Survey 
dFrieman et alJ [2008) toward the end of the 2006 observing 
season at R.A. = 22 h 08 m 53!56, Dec = +00°53'50"4 (J2000). 
Labeled interna lly as a 'strange hostl ess transient', it was ini- 
tially classified (Stritzin ger et alj200 6t) as a. possible SN lb based 
on a spectrum obtai ned at the No r dic Op tical Telescope (NOT). 
In the analysis by Ost man et alJ d201 lh . it was given a SN II 
designation, noti ng, however, that it ma tched less than five tem- 
plates in SNID dBlondin & Tonrvll2007l) . Today, we know that 
these misclassifications were due to the lack of suitable com- 
parison spectra in the literatu re, and, here, we id entify it as an 
H-poor SLSN, as defined bv lOuimbv et al.l d201lh . Our obser- 
vations and results are presented in the next section. Section [3] 
contains the discussion and Section|4]our concluding remarks. 

2. Observations and results 

Photometric observations were carri ed out with the SPSS te le- 
scope at Apache Point Obser vatory dGunn et al.ll 19981 l2006h in 
the SDSS ugriz filters dFukugita et aljl996h . As for all SNe dis- 
covered b y the SPSS SN sur vey, the identification was made ac- 
cording to Sak o et alJ d2008h and the photomet ry was performed 
in the way described in Holt zman et al. (2008),' to which we re- 
fer the reader for more details. We just poin t out th at these are 
'cLsinK magnitudes, defined bv lLupton etalld 1999b to be iden- 
tical to the traditional astronomical magnitudes at higher signal- 
to-noise ratio, but which provide a well-behaved finite flux and 
error even in the low-flux regime. We have derived more tradi- 
tional 3er upper limits (see Table [3]) for non-significant detec- 
tions (error > 0.5 mag), but used the asinh magnitudes for the 
construction of our bolometric light curve (Sect. l3~2b . 

The SPSS light curves are plotted in Fig.Q] As can be seen, 
our observations cover the rise of SN 2006oz, and end before or 
close to maximum light. Because of a six day gap in the obser- 
vations, the explosion date is not very well constrained, but the 
g-band limit previous to our first detection suggests that the ex- 
plosion occurred between MJP 54022 and 54028. Nevertheless, 
these ugriz light curves are very interesting: there are hints in 
our photometry that the SN initially faded in the u-band and re- 
brightened a few days later. This behavior has only been ob- 
served for a handful of stripped core-collapse SNe that were dis- 
covered at very early phases and is usually attributed to a phase 
of adiabat ic cooling of the envelope following sho ck breakout 
(e.g. iRichmond et alj[l994l; IStritzinger et al.ll2002b . Our obser- 
vations are only poorly constraining (because we are in the low- 
flux regime), but based on the error bars and upper limits we 
estimate that there is a probability of 78% that the SN faded 
between MJP 54028 and 54037 in the u-band. Additionally, a 
(significant) dip in the light curve is observed simultaneously 
in the g- and r-bands. Combining the ugr-bands there is strong 
evidence (>5er) for a non-standard, non-increasing luminosity 
evolution of this SN during the first ten days of observations. 

At MJP 54061.87, we obtained a spectrum of SN 2006oz 
at the NOT equipped with ALFOSC. The spectrum spans the 
wavelength range 3200-9100 A with a resolution of 20 A and 
was reduced in the manner described by Ostman et al] d201 lh . 
The spectrum is displayed in Fig. [2] with a sample of compar- 
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Figure 1. ugriz light curves of SN 2006oz obtained by the SDSS SN Survey, shifted as indicated for presentation purposes. 
Triangles denote 3<r upper limits. The lower panel focuses on the g-band evolution shown together with a cubic polynomial fit 
to the last 10 observations (solid line) extrapolated backward (dashed line). Clearly, the first points are strongly inconsistent with a 
smooth rise. The date of the NOT spectrum is also indicated. 



ison spectra from Quimbv et Id] d201 lb . obtained at compara- 
ble pre-maximum phases, to highlight their similarity. To obtain 
a redshift estimate, our spectrum was initially cross-co rrelated 
with the fo ur objects presented in lOuimbv et al.l d201 lb by us- 
ing SNID (Blo ndin & Tonrvll2007l) . The constructed templates 



matched reasonably well, all within 0.365 < z < 0.397 with a 
mean z — 0.376 ± 0.014. It is possible, however, that there are 
intrinsic velocity differences between the features of SN 2006oz 
and the other SNe, which would influence our redshift esti- 
mate. We have searched in the wavelength window ^3500- 
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Figure 2. Spectrum of SN 2006oz ( red) obtained a t the N OT. For 
comparison we show 3 spectra from Ouimbv ~et alJd201 lb at sim- 
ilar pre-maximum phases. The colored spectra are smoothed ver- 
sions (moving average of 5 pixels) of the original gray spectra. 
All spectra were plotted in their rest frame assuming z = 0.376 
for SN 2006oz. The inset shows a zoom of the area around 
2800 A. At this redshift we identified the only probable line (al- 
though at low significance 2.1a) consistent with the Mg II dou- 
blet (which would appear blended at this resolution). 



4500 A for any absorptio n lines, in particular Mg II AA2795, 
2802 (Ou imbv et alj|201 lb that would allow us to determine the 
redshift more accurately. At the low resolution of our spectrum 
the doublet will appear blended, making the search more com- 
plicated. The only candidate line that we were able to identify 
is at ^3856 A. Although at low significance (~2.7<r), it appears 
consistently in all sets of extractions that we attempted and al- 
ways with the same profile (Fig. |2 inset). Although this line is 
insufficient to derive a robust spectroscopic redshift (i.e. Mg II at 
z = 0.377), it can be considered as supporting evidence for the 
redshift derived by cross-correlating the SN spectra. 

Our best redshift estimate, used throughout the paper, is 
therefore z = 0.376. Assuming a cosmology with Hq = 71 km 
&- 1 Mpc" 1 , n m = 0.27 and n A = 0.73, SN 2006oz reached an 
absolute (rest-frame) magnitude of M u — —21.5 mag. This cal- 
culation contains the assumption that maximum light occurred 
in our last epoch of observations and is, therefore, only a lower 
limit to the peak luminosity of the SN. In addition, owing to 
the uncertainty in the cross-correlation redshift, one needs to as- 
sign a 0. 1 mag systematic error in all quantities that depend on 
the distance to the SN. Similarly, we are able to place a strict 
constraint on the (rest-frame) rise time of >29 days. Rest-frame 
wavelengths and related quantities will suffer from an uncer- 
tainty of 1%. 

On 26 August 2011, almost 5 yr after the SN, we obtained 
deep gri imaging of the field with OSIRIS on Gran Telescopio 
Canarias (GTC). The SN position is spatially coincident with a 
candidate host galaxy (object A; Fig. [3]), but within a radius of 7" 
there are four more galaxies, including one that appears extended 
(object E). The photometry of these galaxies and their angular 
distances from the SN position are presented in Table Q] The er- 
ror of these angular distances as determined from the rms of the 




Figure 3. Left: image of SN 2006oz obtained with the SDSS 
telescope (MJD 54061.17). Right: deep i image of the field, ob- 
tained almost 5 yr after the SN with GTC/OSIRIS. The SN is co- 
incident with the position of galaxy A. However, within a radius 
of 7", there are four more objects, including one that appears ex- 
tended, labeled B, C, D and E, respectively. Their photometry is 
reported in TableQ] Both images are 72" x 72". North is up and 
east is to the left. 



Table 1. Host galaxy candidates for SN 2006oz. 



id 


9 


r 


i 


dist. 




(mag) 


(mag) 


(mag) 


(arc sec) 


A 


25.74 ± 0.19 


24.43 ± 0.06 


24.14 ±0.12 


0.47 


B 


> 26.47 


25.62 ±0.16 


25.02 ± 0.25 


2.97 


C 


> 26.47 


25.81 ± 0.19 


24.65 ±0.18 


6.26 


D 


25.18 ± 0.11 


24.70 ± 0.07 


24.32 ±0.12 


4.56 


E 


25.28 ±0.15 


23.71 ± 0.03 


23.21 ± 0.06 


6.56 



Notes. The indexing follows the notation in Fig. [3] All magnitudes are 
given in the AB system. The last column contains the angular distance 
from the SN position (as determined by a geometrical transformation 
between images) and has an associated error of 0.35". 



IRAlQ task geomap was estimated to be 0.35". Because of its 
spatial proximity, we identify galaxy A as the host of SN 2006oz. 
Assuming z — 0.376, the SN occurred at a distance of 2.40 ± 
1.36 kpc from the galaxy center, while the distance to the near- 
est other candidate would exceed 15 kpc. At this redshift, the 
host of SN 2006oz is intrinsically faint with a (rest-frame) ab- 
solute magnitude of M g = —16.9. We note that the photomet- 
ric redshift obtained for this galaxy (u sing the code Le Phare; 
lArnouts et aDll999t Elbert et alJ l2006) is z phot = 0.37to 7 l- 
Although not constraining, this value is consistent with what was 
obtained for the SN. 

3. Discussion 

3.1. Black-body fits 

To estimate the photospheric temperature and radius of 
SN 2006oz, we fitted a black body (B B) function to the op - 
tical spectrum. As also pointed out by Chomi uk et alJ d201 lb . 
these events miss some flux in the UV with respect to a BB, 
and therefore reliable BB fits can only be obtained by fitting 
data redward of rest-frame 3000 A (i.e. the Rayleigh-Jeans tail). 
This was also the experience we gained by fitting the NOT 
spectrum. We achieved a reasonable fit with Xdof = 2. 08 over 
1309 data points, resulting in T BB = 15540 ± 430 K and 

1 IRAF is distributed by the National Optical Astronomy 
Observatory: |http : //iraf ■ noao . edu/ iraf /web/1 
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Table 2. Black-body fits to SN 2006oz. 



MJD 


Tbb 


Rbb 




(K) 


(10 15 cm) 


54028.17 


13960 ±460 


0.94 ± 0.05 


54030.18 


15660 ± 3780 


0.93 ± 0.30 


54032.17 






54035.25 


9253 ± 2775 


1.73 ± 1.04 


54037.16 


1 1020 ± 3594 


1.44 ±0.83 


54040.17 


22610 ± 15485 


0.71 ±0.50 


54047.16 


11230 ± 2425 


2.62 ± 0.92 


54050.16 


14420 ± 350 


2.00 ± 0.07 


54053.16 


16270 ± 770 


1.88 ± 0.12 


54056.16 


19160 ± 3280 


1.61 ±0.32 


54057.13 


15890 ± 1170 


2.13 ±0.21 


54059.16 


14820 ± 540 


2.44 ±0.12 


54061.17 


13740 ± 1030 


2.81 ±0.30 


54061.87" 


15540 ± 430 


2.18 ±0.08 


54062.16 


16470 ± 970 


2.19 ±0.16 


54068.15 


15640 ± 3940 


2.50 ± 0.80 



Notes. ' a ' This epoch is a fit to the spectrum. 



Rbb = (2.18 ± 0.08) x 10 15 cm. It is also possible to fit the 
broad-band photometric data, although with higher uncertain- 
ties. We therefore performed a weighted fit of the griz data, ex- 
cluding the tt-band that corresponds to rest-frame 2574 A. The 
fit results are shown in Tableland Fig. [4] and indicate that the 
emitting radius increased from ~1 to > 2 x 10 15 cm during 
the period of our observations. A linear increase describes the 
Rbb evolution well with a weighted fit resulting in Xdof = 1-21. 
In contrast, the temperature does not show any significant evolu- 
tion: a linear fit results in a shallow slope, while fitting a constant 
temperature (—15000 K) is also consistent with the observations 
(Xdof ~ 1-92 versus 1.27). However, during the first ten days the 
average temperature might have been lower (—12400 K). The or- 
ders of magnitude obtained and the Rbb evolution are consistent 
with what has been reported p reviously, although a cooling of the 
BB is usually observed (e.g. Ouimbv e t alj|201 ltlChomiuk et alj 
1201 lb . However, cooling is only observed conclusively at post- 
maximum phases, while, as discussed below, we are probing ear- 
lier stages of SLSN evolution with SN 2006oz. 
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Figure 4. Black-body fits to the SN 2006oz data. In the top 
panel, a fit to the NOT spectrum is shown. Only data red- 
ward of 3000 A (rest-frame) are included in the fit. The middle 
and bottom panel show the radius and temperature evolution of 
SN 2006oz assuming it can be described by a black body. The 
filled circles are the results of fitting the griz photometry, while 
the spectrum fit is represented by a square symbol. The dashed 
lines show linear weighted fits to the data. 



3.2. Bolometric and absolute light curves 

To construct a (pseudo-) bolometric light curve of SN 2006oz 
we integrated the flux included within the limits of our simul- 
taneous ugriz observations. To account for the epochs con- 
taining upper limits (in u and z) we use d the asinh magni- 
tudes provided by th e SDSS photometry dLupton et al.l 1999: 
iHoltzman etaT]|2008l) . e.g. 2.6±1.5 ^Jy for our u observation 
at MJD 54037. This is the most reliable assumption and the high 
uncertainties are propagated to the final bolometric light curve. 
A correction was appli ed only for foregroun d Galactic extinction 
(E{B -V) = 0.047: ISchlegel et al.lll998h but not for any host 
galaxy extinction. The bolometric light curve is plotted in Fig. [5] 
(upper panel) together with that of PSl-lOawh (z = 0.9084), 
which is the only other H-poor SLSN with a pre-maximum bolo- 
metric light curve (based on multi-filter observations) extending 
back to more than a week before maximum light. We constructed 
the bolometric light curve ofPSl-lOawh followin g the same pro- 
cedure as for SN 2006oz based on the data of IChomiuk et alj 
d201 lb . but did not apply a bolometric correction for the red 
tail of the BB (this is why our light curve might appear a lit- 



tle different from theirs). Determining a bolometric correction 
is not trivial. The rest-frame wavelengths probed were —2360- 
6980 A (and -2270-4720 A for PS 1 - 1 Oawh; extending the inte- 
grations out to 0.5 x FWHM from the central wavelengths of the 
bracketing filters). Assuming the limiting case that these SLSNe 
can indeed be described by a BB, we can roughly estimate the 
bolometric correction by measuring the area below the BB curve 
covered by our observations. For the BB fit to the spectrum of 
SN 2006oz, the probed wavelength window captures 51% of 
the emitted flux, while another 43% is at bluer wavelengths and 
only 6% in the red tail. The corresponding percentages for PSl- 
lOawh are 40%, 44%, and 16%, respectively. These rough num- 
bers illustrate that the bolometric luminosities obtained could be 
subject to a correction of up to a factor of 2-2.5. This is the worst 
case, however, because the flux in th e UV clearly deviates from 
a BB (see also lChomiuk et al.ll201 II) . 

In addition to the uncertainty of the total bolometric lumi- 
nosity, there is also an uncertainty of the time of maximum. This 
is illustrated in Fig. [5] by applying a horizontal shift to the LC 
with respect to the limiting case for which we assume that our 
last observation was exactly at maximum light (vivid and faded 
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of PSl-lOawh, which correspond to the same rest-frame wave- 
lengths. This comparison shows that (despite the uncertainty in 
tmax) the SDSS observations of SN 2006oz reach deeper limits 
in an absolute level and provide a window to the early evolution 
of an SLSN. 

The shape of this very early bolometric light curve is very 
interesting. The individual band observations discussed in the 
previous section were translated into a plateau-like phase in the 
bolometric light curve. This plateau has a minimum duration of 
six days and is followed by a dip, after which the luminosity be- 
gins to rise. Owing to the uncertainty in the time of explosion, 
the duration of this initial phase can be up to 14 days. Such a 
precursor plateau has not been observed before for an SLSN, so 
the reasonable question is whether it is intrinsic to all similar ex- 
plosions or special to SN 2006oz. Its detection was possible by a 
combination of the SN redshift with multi-color observations in- 
cluding blue filters, a criterion not met by other H-poor SLSNe. 
The p re-maximum observations of the PTF SNe (Ouimbv et al. 
2011) were mostly obtained in the i?-band and d o not extend to 
the res t-frame UV, while the PS 1 observations of lChomiuk et al.l 
(1201 lh are not as deep, due to the longer distances (Fig. [5j. For 
SCP06F6 (Ba rbarvetalJl2009l) . there might be hints for such a 
plateau in the z-band but the photometric errors are too large to 
allow any firm conclusion. It is therefore an intriguing possibil- 
ity that this behavior is intrinsic to all events in this class. It is 
unlikely that this result will be changed because this bolometric 
light curve by construction misses some flux in the UV. Instead, 
the result will more probably be enhanced because the colors of 
SN 2006oz evolve from blue to red during this period, which re- 
duces the fraction of missing flux with time. The nature of the 
plateau is discussed in more detail in Sect. 13.41 



Figure 5. Upper panel: Bolometric light curve of SN 2006oz 
(circles) compared to PSl-lOawh (solid line; data from 
IChomiuk et al.l I20TTI) . which has similar pre-maximum multi- 
color observations. To illustrate the uncertainty in t max of 
SN 2006oz, its light curve was shifted horizontally with respect 
to the limiting case, where our last observation was assumed to 
be taken exactly at maximum light. Assuming the two events 
are similar, we obtain a reasonable match of the light curves by 
shifting t max by —10 days (faded symbols). Middle and lower 
panels: The absolute pre-maximum light curves. The absolute 
luminosities at rest-frame ^4530 A (observed r- and z-bands) 
and ^2560 A (observed u- and g-bands) are compared, after cor- 
recting for cosmological expansion. As above, vivid and faded 
colors show the light curves of SN 2006oz shifted by 10 days. 
A few (non-constraining) upper limits were removed for presen- 
tation purposes. The early observations of SN 2006oz are below 
the detection limits of the PanSTARRSl observations of PSl- 
lOawh. 



colors). This experiment was made to show that it cannot be es- 
tablished whether SN 2006oz was a fainter event or whether we 
are probing earlier stages of its evolution (or a combination of 
both). Indeed, a reasonable match of the light curves can be ob- 
tained by shifting t max by —10 days (rest-frame, time-dilation 
corrected). It is certainly true that the very early data-points of 
SN 2006oz were below the sensitivity of the PanSTARRSl ob- 
servations of PSl-lOawh at z ~ 0.9. This can also be seen in 
the lower panel of Fig. [5] which contains a comparison of ab- 
solute rest-frame light curves. In particular, the observed u and 
r -bands of SN 2006oz are compared to the g and z observations 



3.3. The spectrum 

The lack of clear hydrogen features is an import ant clue to the 
natur e of this event and the others in its class (Ouimbv et al. 
1201 ll) . It is difficult to avoid the conclusion that the object is 
hydrogen-deficient, although the possibility that some hydrogen 
is present, but difficult to detect, cannot be rejected. The lack 
of clear signatures of helium is a more ambiguous clue. The 
prominent lines of He I in the optical and NIR arise in tran- 
sitions to the n = 2 level from higher levels. Ambient con- 
ditions in SN photospheres are typically insufficient to popu- 
late these upper levels, consequently they must be populated 
by non-thermal dHarkness etalJll987r. iLucvll 19911) or non-LTE 
dDessart et al.l l201 1 ) processes. Relevant ionization energies are 
H - 13.6 eV, He I - 24.6 eV, C I - 11.3 eV andOl- 13.6 eV. From 
the absence of obvious helium lines, it is then much less obvious 
that an absence of helium can be firmly deduced. In addition, the 
wavelength range of the SN 2006oz spectrum (and most similar 
events), does not cover the most prominent He lines, with the ex- 
ception of A5876 , which is found in a noisy part of the spectrum. 
The question whether He might be expected requires additional 
investigation and better wavelength and temporal coverage. 

We have app l ied the parametr ized code SYNOW 
(lHatano et al.l 119991: iBranch et afl 120031) to model the spec- 
trum of SN 2006oz. Although SYNOW is based on simplified 
physical assumptions, it is very useful for identifying spectral 
features that are strongly Doppler-broadened and sometimes 
heavily blended. Because SYNOW contains many adjustable 
parameters, we did not attempt an automated fitting via chi- 
squared minimization but instead looked for a reasonable 
agreement between the observed and the model spectra by eye, 
using the minimum number of different atoms/ions in the enve- 



6 



:r-luminous SN 2006oz 



10 ■ 



SCP06F6 
SN 2006OZ 
SYNOW w c n 
SYNOW w/o C II 



fa 
13 

■a 



i 

•a 



o 



0.1 




2000 



5000 



6000 



3000 4000 _ 

Rest Frame Wavelength (A) 

Figure 6. SYNOW fit to the spectra of SN 2006oz and SCP06F6. 
The observed spectra were brought to the rest-frame by assum- 
ing z = 0.376 and z = 1.189, respectively, and their fluxes 
were scaled to match in the overlapping wavelength region. Two 
SYNOW models, with and without C II, are overplotted. The inset 
shows a zoom in the spectrum of SN 2006oz around 3300 A. The 
magenta line (inset only) shows a model without Fe III that pro- 
vides a worse fit in this region. The lower panel shows the con- 
tribution of the individual ions to the model spectra. The best- 
fitting model has r(0 II) = 0.02, r(Mg II) = 0.05, r(Si III) = 0.5, 
r(S III) = 0.5 and r(Fe III) = 0.4. 



lope. The photospheric expansion velocity and the temperature 
of the underlying blackbody radiation were initially set and kept 
fixed during the search for features. Adopting v p \ iot = 12000 
km s _1 and Tbb = 14000 K resulted in a good fit to most 
spectral features. A power-law atmosphere, where the optical 
depth as a function of velocity varies as {v/v max )~ n with 
n = 7, has been assumed for all atoms/ions. Our tests showed 
that there are no significant differences between the line profiles 
obtained by using either power-law or exponential optical depth 
profiles for the spectral features of SN 2006oz. The maximum 
velocity of the envelope u max was set to 40000 km s _1 , but this 
parameter is only weakly constrained by our spectrum. 

The observed spectrum and two model spectra, with and 
without C II, are plotted togeth er in Fig. [6] The figu re also con- 
tains the spectrum of SCP06F6 (iBarbarv et al.l2009t) to illustrate 
a possible extension of the spectrum to the (rest-frame) UV. The 
contributions of the individual ions to the model spectrum are 
also shown. The "W"-like feature aroun d 4200 A nicely match es 
the O II feature that was also invoked bv lOuimbv et al.l d201 ll) to 
model other SLSNe. The strong feature around 2700 A is prob- 
ably due to Mg II, although the spectrum is noisy in that region. 
Fein can explain the observed bump around 3180 A (Fig. |6j 
inset), and none of the other ions considered here were able 
to model that feature. Si III has a contribution around 3000 A, 



but its presence was mostly motivated by the assumption that 
SN 2006oz is spectrally similar to other SLSNe. Indeed, the 
strongest Si III feature matches the obser ved feature aroun d 
2420 A in SCP06F6 and in other SLSNe (lOuimbv et al.ll201ll) . 
We have also tried to include S III in the model of SN 2006oz, but 
its spectral features are weak and easily obscured by the noise. 
We consider the presence of Si III and S III as a possibility, but 
they cannot be directly identified in SN 2006oz. 

Th is model is consistent with the findings by lOuimbv et al.l 
d2011h concerning the presence of O II, Mg II and Si III (assum- 
ing there is indeed some similarity among the spectra of these 
extreme SNe). In addition, we suggest that Fe III might also be 
present and it may also be a good candidate to model the sharp 
drop of the observed flux around 2000 A in SCP06F6. We did 
not find the presence of C II necessary to model SN 2006oz. This 
may be partially because the strongest f eature that was explained 
with this ion by Ouimb v~et alJ d201 ll) is beyond the observed 
spectral range (rest-frame 2200 A). On the other hand, Fig. [6] 
illustrates that C II has strong contributions around 2800, 3800, 
and 4500 A, but these features do not match the observed spec- 
trum. There is some ambiguity between the presence of C II and 
that of O II plus Mg II, but the latter ions better match the ob- 
served features than C II alone. We also note that for C II, unlike 
for all other ions for which T oxc = Tbb was assumed, we had to 
apply a lower excitation temperature (10000 K), because hotter 
temperatures produced unacceptably strong features. 

The chemical composition identified above (O II, Mg II, and 
maybe Fe III, Si III and S III) is consiste nt with the "carbon - 
burned" SN atmosphere considered by Hata no et al.l (l999). 
Although it is probably premature to conclude that SN 2006oz 
had such an atmosphere, it is interesting that together with the 
lack of H and He I and maybe C II in the spectrum, all identified 
features belong to the ions that are expected to be the strongest 
at an excitation temperature of ^14000 K in an atmosphere con- 
taining elements that underwent carbon-burning. 

3.4. Models for the light curve 

The bolometric light curve can be divided into two parts: the 
initial precursor plateau and the subsequent smooth monotonic 
rise. As discussed above, it is not clear if the maximum of the 
light curve was observed, a handicap in constraining models. 

We note that precursor pla teaus, some similar i n shape to that 
of SN 2006oz, were found bv lDessart et al.l d201 ll) in their mod- 
els of helium star explosions. In their models, the plateau is as- 
sociated with the shock breakout, fireball-cooling phase, before 
the rise to the peak powered by radioactive decay. The plateau 
arises when the temperature decrease at the photosphere slows 
down, an effect associated with the recombination of ejecta lay- 
ers to their neutral state (primarily He, but also CNO elements 
in the helium-rich progenitor models). The plateau brightness is 
determined by the amount of energy initially deposited by the 
shock and the size of the progenit or envelope. The bolometric 
light curves of Des sart et alJ (1201 ll) typically have a plateau that 
lasts about ten days and is about a factor of 20 to 30 dimmer than 
the subsequent peak. 

The physics e xplored by the precursor plateaus of 
iDessart et alJ d201 ll) may be relevant to SN 2006oz and its kin, 
but the plateau these authors find seems to be too dim to directly 
correspond to that observed in SN 2006oz. In SN 2006oz, this 
plateau is dimmer than the subsequent peak by about a factor 
of 8-10 in luminosity (depending on £ max ; Fig. [5]). Becau se th e 
plateau is brighter than in the models of Dessart et aj] d201 lb . 
the material must be distributed at larger radii. The structures 
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Figure 7. Bolometric light curve of SN 2006oz shown along 
with models describing the rise to maximum (ignoring the pre- 
cursor plateau). Models include the following power sources for 
the luminosity: radioactive decay (red, dashed), input from a 
magnetar (green, dot-dashed) and CSM interaction (blue, solid). 
Details concerning the model parameters are given in the text. 



in the models of Dessart et aT] (1201 ll) are originally in hydro- 
static equilibrium as dictated by the systematics of stellar evo- 
lution. A substantially larger radius suggests that the plateau we 
observe in SN 2006oz does not arise in a stellar envelope, but 
in a circumstellar medium. Consequently, the only likely expla- 
nation that we have found for the plateau is that it might rep- 
resent a recombination wave in a CSM that surrounds the pro- 
genitor star. A possible case of shock b reakout from a H-rich 
dense CSM was discussed bv lOfek et afl(l2010b . For SN 2006oz, 
unlike for SNe IIP, this plateau cannot be caused by recombi- 
nation of H, which occurs at much lower temperatures (5000- 
6000 K). The derived temperature also makes a He recombi- 
nation explanation unlikely, because most of the He atoms are 
singly ionized at T > 10000 K, but dou ble ionization starts 
above T > 15000 K (lHatano et al.lll999l) . On the other hand, 
a recombinati on plateau could be consistent with the transition 
Olll to Oil (lHatano et all 1 19991) . This idea is also potentially 
consistent with the detection of O II in the spectrum of H-poor 
SLSNe (Sect. I3.31 >. Interesti ngly, a precursor plateau was ob- 
served in the simulations of Blinnikov & Soro kinal (1201 Ol) for 
SN la explosions in a C-0 CSM, but, as the authors argue, it was 
an artifact of the arbitrary initial conditions. An early plateau 
is also pres ent in the bolome tric light curve of the 1987A-like 
SN 2006au dTaddia et al.ll2012l) . 

The subsequent smooth rise of the light curve represents a 
more normal SN behavior and could be satisfied by a variety of 
processes, such as radioactive decay, input from a pulsar or mag- 
netar, or luminosity from SN ejecta-CSM interaction. General 
se mi-analytic models of these processes have been constructed 
by IChatzonoul os et all (120121) . and will be applied to a range 
of superluminous events in Chatzopoulos, Wheeler & Vinko (in 
prep.). Figure [7] shows a sample of models compared with the 
observed bolometric light curve of SN 2006oz, ignoring the pre- 
cursor plateau. The models all peak at a time around the last 
observed data point, but the lack of a measured maximum (and 
subsequent decline) means the models are not well constrained. 



The models were not fitted in any formal sense, but constructed 
to be reasonable "by eye." The ejecta mas ses derived may be 
somew hat underestimated as discussed by IChatzopoulos et all 
(12012b . 

The red dashed curve shows a model driven by radioactive 
decay with M^i = 10.8 Mq and a diffusion time = 36 days, 
determined by the time scale of the rise of the light curve. The 
SN ejecta mass, M C j = 14.4 Mq, is derived from the diffusion 
time, ta, assuming a velocity, v = 12000 km s _1 , and opacity 
k = 0.05 cm 2 g -1 . The ejecta mass scales as k^ 1 and conse- 
quently would be lower for higher opacity values. The fact that 
such a big fraction of Af c j has been burned to 56 Ni, unlike in 
most observed SN types, shows that a radioactive decay model 
has difficulties to account for the rise of SN 2006oz to maximum. 
Furthermore, the SN was not detected in the next observing sea- 
son that started nine months later. By stacking the images from 
the first five weeks, we deduce a limit r > 24.2 (Table |3]l im- 
plying a decay rate > 2.03 mag 1Q0 -I d — ^ This is steeper than 
the radioactive decay of 56 Co (0.98 mag 100 _1 d _1 ), assuming 
complete 7-ray and positron trapping. The radioactive model for 
H-poor SLSNe has also been ser i ously challenged by th e ob- 
servations of iQuimbv et al.l (1201 lb . iPastorello et alJ d2010l) . and 
IChomiuk et al.l (1201 lb . We stress, though, that an assumption in- 
volved i n these calc ulations is that that 56 Ni is centrally con- 
densed dArnettl[l982b . It is therefore possible that a configura- 
tion where 56 Ni is substantially mixed into the outer envelope, 
thereby allowing 7-rays to escape freely, cannot be ruled out. 

For the same velocity, opacity, and rise time, a magnetar 
model yields the same ejecta mass, 14.4 Mq, as the radioac- 
tive decay model, but requires two more parameters, E p and t p , 
where the decay model requires only one, Msa, to fit the rise and 
peak. The green dot-dashed curve shows a magnetar spin-down 
model with initial rotational energy E p — 1.45 x 10 51 erg and 
a characteristic spin-down time of t, p = 13 days, giving an initial 
rotational period 3.7 ms and a magnetic field B = 2.24 x 10 14 G. 

The blue curve shows a model with SN ejecta and CSM in- 
teraction (including both forward and reverse shocks) for again 
the same velocity, opacity, rise time and ejecta mass as the pre- 
vious models. This CSM model was adjusted to give both a 
reasonable fit to the rise time and a reasonable formal black- 
body temperature at maximum . The other parameters used are 
("see IChatzonoulos et al.ll20l2l for detailed definitions): E sn — 
2.5 x 10 51 erg, the density slope of the SN ejecta n = 12, 
the progenitor radius R p = 2.5 x 10 14 cm, the density slope 
of the CSM s — 0, the density scale at the base of the CSM 
Pcsm.i = 5 x 10~ 13 g cm~ 3 (determined at R p and equal to 
the constant density of the CSM shell), M csm = 6.5 Mq, and 
MNi = 0.02 Mq. The constant CSM photospheric radius of this 
model is i? p h = 1.8 x 10 15 cm, virtually the same as the radius of 
the shell itself, giving an optical depth of r csm = 327. For these 
parameters and for L^ \ — 1.3 x 10 44 erg s _1 at peak, the formal 
black-body temperature at maximum light is Tbb = 15300 K. 
This model, with several free parameters, can also provide a de- 
cent fit to the data. 

The uncertainty in the bolometric correction and t max has, 
of course, an impact on the models discussed above, although 
this is not easy to quantify because of parameter degeneracies. It 
is, however, almost certain that constructing a viable radioactive 
model will become increasingly difficult as these uncertainties 
can only increase the necessary M-^i (and the ratio Mjsn to M e j) 
to power the light curve. Our experience shows that it will still be 
possible to obtain viable magnetar or CSM models by modifying 
parameters such as E p , t p , p csm and M csm . More constraining 
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for these models would be the availability of post-maximum data 
and decay rate. 

We have shown that a range of models can explain the 
smooth rise to maximum. Combining these models with the pre- 
cursor plateau is less straightforward, as none of these models 
are presently able to numerically reproduce it. The possibility 
that the plateau is the result of a recombination wave in a (O- 
rich?) CSM, the only reasonable suggestion we have come up 
with, gives a rationale for considering that the same shell might 
be responsible for the rise. This picture is not perfectly self- 
consistent, however, because the first part requires the shock to 
have broken out of the CSM, while the second part requires the 
shock to still be interacting with it. It is therefore difficult to 
understand how the CSM was heated in the first place, i.e. be- 
fore the ejecta-CSM interaction provides the rise to the peak. A 
possible answer could be by energy deposition to the CSM by 
the SN blast wave, although this possibility remains to be stud- 
ied. Another alternative solution could be provided by a hybrid 
model with a CSM recombination accounting for the plateau and 
a magnetar for the rise to the peak. Additional observations of 
similar events, including more multi-color light curves from very 
early to very late phases, will be required to resolve this problem. 

3.5. Host Galaxy 

The host of SN 2006oz is only the fourth out of nine sim- 
ilar H-poor SLSNe to have been identified dNeill et al.ll201lt 
IChomiuk et al.|[20TTh . Very little information is available on the 
hosts of these events, and many of them are detected in only one 
filter, so that a study of the properties of the host of SN 2006oz 
is warranted. 

The colors g-r = 1.26 ± 0.20 and r - i = 0.26 ± 0.13 
(after correcting for Galactic foreground extinction), suggest a 
significant break between the g and r filters. This can only be as- 
sociated with the 4000 A break, at a redshift consistent with the 
one we are examining. These colors are indeed consistent with 
an intermediate value D n (4000) = 1.4 ± 0.1 for the 4000 A 
break, indicative of relatively old stellar populations and, pos- 
sibly, of a burst of star formation within th e last 1-2 Gyr (e.g. 
iKauffmann etai]|2003t iGallazzi et alj |2005|). In Fig. [8j the host 
galaxy colors (if -corrected to z — O.ljjare compared to a set of 
galax ies, selectecjl from the main spec troscopic sample of SDSS 
DR4 dAdelman-McCarfhv et al]|2006l) to have 0.08 < z < 0.12, 
and to a set of du st-free model galaxies w ith different star- 
formation histories dBruzual & Ch ariot 2003). Despite the high 
uncertainty, the galaxy position on this color-color diagram is 
more consistent with low-mass, star-forming galaxies and is far 
from the locus of elliptical, passive galaxies. It is not, however, 
a blue , starburst galaxy: by using the temp late of Kinney et alj 
(1996) and assuming the extinction law of Calzetti et al. (2000), 
it is not possible to reproduce these colors for any value of red- 
dening. 

To gain a more quantitative understanding of the nature of 
the host galaxy, its observer-frame colors were compared to 
those predicted by a Monte Carlo library of model spectra , 
redshifted to z = 0.376, based on iBruzual & Chariot! d2003l) 
population synthesis models spanning a wide range of star 
formation histories, metallicities, and dust attenuation (mod- 



2 The /^-corrections were estimated from the model that best fits the 
observer-frame galaxy colors, as described below, and have an uncer- 
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Figure 8. Color-color diagram comparing the host of 
SN 2006oz (black cross) to SDSS galaxies (blue) and to a set 
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eled as in ICharlot & FalU 120001) . In this way, we constructed 
probability density functions and derived the following esti- 
mates for the host properties: a stellar mass of log(M /Mq) = 

8.7 ± 0.23, a luminosity-weighted age of 1.28^0.61 G Y r ' and a 
SFR of 0.17l„ J* Mq yr" 1 . This results in a specific SFR of 
log(SFR/M*) = -9.47 ± 0.39 yr " 1 , consistent with , but at the 
lower end of, the values reported bv lNeill et al.1 d201 ll) (although 
their uncertainties are high). 

With an absolute luminosity of M g = —16.9, this galaxy 
is probab ly metal-poo r. By using the mass-metallicity rela- 
tion of iTremonti et al.l (2004), we obtain an indicative value 
of 12 + log (O/H) = 8.34 ± 0.16 for the oxygen abundance 
th at corresponds to 0.4 5 Zq (assuming the solar abundances 
of lAsplund et al.l [20091) . A direct met allicity measure ment ex- 
ists only for the host of SN 201Qgx dStoll et al.ll201 ll). indicat - 
ing 0.47 Zq dKobulnickv & Kewlevl2004t lAsplund et al.l2 009). 
These values are not particularly low, but the fact that all H-poor 
SLSNe have been found in faint galaxies suggests that metal- 
licity might indeed play a role. The reason why this is interest- 
ing is that all viable models leading to H-poor SLSNe involve 
(one way or another) a very massive star that has suffered ex- 
tensive mass loss. Since stellar winds ar e driven by metals (e.g. 
IVink & de Koterll2005l: iPuls et alj|2008h and are less efficient at 
low metallicity, it has been suggested that SLSN progenitors lose 
mass through episodic outbursts due to mechanisms that are not 
well understood, possibly involving pulsational pair instability 
(e.g.[Wooslev et al. 2007). On the other hand, a similar discus- 
sion has been going on about the preference of long GRBs, with 
optical afterglows, to be found in metal-poor galaxies, but it has 
been shown that this might b e partly due to the fact that sm aller 
galaxies have a higher SFR (Ma nnucci et al.l I2010L 120111) . It is 
therefore possible that the same bias exi sts for the SLSN hosts . 
The fundamental metallicity relation of iMannucci et al .1 (1201 lh 
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for low-mass galaxies predicts 12 + log (O/H) = 8.39 for the 
host of SN 2006oz, i.e. again a value that is not very low. 

4. Conclusions 

We have studied SN 2006oz , an event belonging t o the family of 
H-poor SLSNe identified bv lOuimbv et alJ (1201 lh . The bolomet- 
ric light curve shows a precursor plateau with a rest-frame dura- 
tion of 6-10 days, while the bluest bands demonstrate a dip in the 
luminosity, before rising smoothly to maximum light. Our early 
observations are more sensitive than those of previous studies 
and it is therefore possible that this behavior might be common 
to more H-poor SLSN. If this is the case, it provides an impor- 
tant diagnostic for the nature of these events. We argued that the 
precursor plateau might be caused by a recombination wave in 
an H-deficient CSM. The subsequent rise in luminosity can be 
described by energy input from a magnetar or by ejecta-CSM 
interaction, but radioactive decay is less likely. Although the 
model parameters are not well constrained owing to the lack of 
post-maximum data, all models involve large ejected and CSM 
masses (if a CSM is required), pointing to a very massive pro- 
genitor. 

Our deep observations with GTC have revealed the faint host 
galaxy of this event (M g = —16.9). It is a moderately low-mass 
(log(M/M ) = 8.7), star-forming (0.17 Af yr _1 ) galaxy with 
a luminosity-weighted age between 0.7-2.3 Gyr. It is not a blue, 
starburst galaxy and the metallicity, inferred indirectly, is not 
particularly low. 
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G. Leloudas et al.: The super-luminous SN 2006oz 
Table 3. SDSS photometry of SN 2006oz (AB magnitudes). 



MJD 


u 


9 


r 


i 


z 


54010.19 


>21.752 


> 24. 124 


>22.759 


>22.539 


>21.128 


54012.18 


>21.803 


>21.650 


>22.360 


>21.805 


>21.321 


54020.18 


>22.740 


>24.200 


>22.673 


>22.688 


>20.332 


54022.23 


>22.228 


>23.674 


>22.898 


>22.772 


>20.815 


54028.17 


21.391 (0.198) 


22.280 (0.112) 


22.314(0.237) 


22.359 (0.388) 


> 20.666 












22.092 (0.856) 


54030.18 


21.823 (0.183) 


21.940 (0.067) 


21.926 (0.116) 


22.268 (0.230) 


>21.293 












22.683 (0.689) 


54032.17 




21.830 (0.217) 








54035.25 


>20.899 


22.591 (0.266) 


22.451 (0.434) 


22.251 (0.425) 


>20.184 




21.763 (0.440) 








21.152 (0.514) 


54037.16 


>21.795 


22.240 (0.141) 


22.238 (0.320) 


21.858 (0.227) 


>21.258 




22.861 (0.595) 








23.096 (0.833) 


54040.17 


>21.050 


21.552 (0.216) 


21.567 (0.232) 


22.299 (0.515) 


> 20.602 




22.490 (0.990) 








22.071 (0.904) 


54047.16 


>20.364 


20.924 (0.219) 


20.596 (0.174) 


20.872 (0.304) 


> 19.567 




21.547 (0.713) 








21.306(1.372) 


54050.16 


20.723 (0.114) 


20.534 (0.032) 


20.590 (0.051) 


20.666 (0.089) 


20.749 (0.378) 


54053.16 


20.337 (0.082) 


20.283 (0.025) 


20.412 (0.040) 


20.507 (0.065) 


20.843 (0.318) 


54056.16 


20.408 (0.172) 


20.160 (0.034) 


20.302 (0.059) 


20.694 (0.145) 


20.571 (0.368) 


54057.13 


20.332 (0.067) 


20.089 (0.021) 


20.150 (0.033) 


20.370 (0.056) 


20.642 (0.225) 


54059.16 


20.097 (0.079) 


20.011 (0.021) 


20.067 (0.035) 


20.177 (0.050) 


20.520 (0.225) 


54061.17 


20.115 (0.079) 


19.964 (0.026) 


19.930 (0.034) 


20.122 (0.052) 


20.098 (0.184) 


54062.16 


20.318 (0.092) 


19.927 (0.030) 


20.017 (0.038) 


20.226 (0.062) 


20.467 (0.285) 


54068.15 




19.800 (0.079) 


19.810 (0.076) 


20.217 (0.161) 


19.909 (0.242) 


54367.02" 
54395.57 6 


>23.228 


>24.508 


>24.209 


>23.592 


>22.190 



Notes. The quoted limits are 3a. For the epochs after th e SN discovery, where the SN was below the formal detection limit (it and z -band), we 
also provide, below the 3er limits, the asinh magnitudes (Lupton et al. 1999; Holtzman et al. 2008) , which were used as the best approximation 
in the construction of the bolometric light curve. Observation mid-point and magnitude limits from stacking all observations from the first 5 
weeks of the 2007 season. ^' Observation mid-point and magnitude limits from stacking all observations from the 2007 season. 
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